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Abstract 


Young women with disabilities face unique barriers in the transition from school to adulthood, 
yet very few studies have examined the effectiveness of gender-specific career interventions. 
Using an intent-to-treat analysis, this study tested the efficacy of the Paths 2 the Future (P2F) 
career development curriculum to produce beneficial impacts as compared to business-as-usual 
career and transition services within a clustered, school-randomized trial. The sample included 
366 young women with disabilities enrolled in 26 high schools. Controlling for student and school 
characteristics, multilevel growth models showed that young women in the P2F intervention 
schools grew in career development skills at a greater rate relative to those in the control 
schools. Findings suggest that P2F was effective in promoting individual student differences in 
career development skills for young women with disabilities and support the need for further 
research examining effectiveness of gender-specific career interventions in special education. 


Although the Individuals With Disabilities 
Education Act (2006) states that the funda- 
mental purpose of a free appropriate public 
education is to prepare all students with dis- 
abilities for “further education, employment, 
and independent living” (20 U.S.C. § 1400[d] 
[1][A]), youth with disabilities consistently 
experience lower rates of employment and are 
less likely to enroll in or complete postsecond- 
ary education than their nondisabled peers 
(Butterworth & Migliore, 2015; U.S. Bureau 
of Labor Statistics, 2020). Young women with 
disabilities, in particular, face unique chal- 
lenges in accessing postsecondary education, 
gaining employment, living independently, 
and fully participating in their communities 
(Hogansen et al., 2008; Lindsay et al., 2017). 
After leaving high school, females with dis- 
abilities are more likely than males to be 
employed part-time and earn lower wages and 
are less likely to work in high-skill or high- 
wage jobs (Doren et al., 2011). The persistence 
of these inequities over time suggests that 


women with disabilities are more likely to 
experience poverty in adulthood (Fins, 2019). 
The process of career development for young 
women with disabilities is impacted by both indi- 
vidual and structural barriers. Individual attributes, 
such as low self-esteem, limited self-efficacy, and 
a lack of self-advocacy skills can restrict the 
ability to fully explore a wide range of career 
options. In addition, young women with disabili- 
ties often have limited opportunities for career 
exploration activities during high school (Ferri & 
Conner, 2010; Lindstrom et al., 2012), are less 
likely than their male peers to enroll in vocational 
courses or participate in community work experi- 
ences, and may face discrimination related to both 
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gender and disability identities (Lindsey et al., 
2017; Noonan et al., 2004). This restricted set of 
opportunities and experiences translates to a nar- 
row range of career interests and aspirations, 
ultimately resulting in poor long-term educa- 
tional and employment outcomes (Doren et al., 
2011; Lindsey et al., 2017). 

Despite the complex combination of barri- 
ers faced by young women with disabilities as 
they transition from school to adult roles, rela- 
tively few career development interventions 
have been developed to specifically focus on 
ameliorating these gender gaps. One of the few 
exceptions is Career Connections, a 24-lesson 
curriculum designed to improve prosocial and 
career development outcomes for young 
women with disabilities. This weekly program 
aims to improve participants’ knowledge about 
self-advocacy, gender equity, leadership, and 
workforce issues and provides opportunities to 
investigate a variety of career options. A pre- 
test-posttest study showed that immediately 
after participating in Career Connections, 
female students with disabilities who partici- 
pated (N = 102) had a range of improved out- 
comes, including (a) more positive relationships 
with adults, (b) increased ability to advocate for 
themselves in difficult situations, (c) increased 
participation in school and community activi- 
ties, and (d) more proactive communication 
skills (Lindstrom et al., 2008). 


Despite the complex combination of 
barriers faced by young women 
with disabilities as they transition 
from school to adult roles, 
relatively few career development 
interventions have been developed 
to specifically focus on 
ameliorating these gender gaps. 


Another such intervention, Girls at Work, 
developed by Wehmeyer and colleagues (2009), 
utilized a self-determined career development 
model to improve career outcomes for young 
women with intellectual and developmental 
disabilities. Girls at Work is a web-based, stu- 
dent-driven, transition-oriented curriculum that 
includes a set of eight action steps designed to 
promote self-determination and encourage gen- 


der awareness. Initial field test data indicated 
improved employment rates and enrollment in 
postsecondary education for the young women 
(N = 54) who completed the Girls at Work pro- 
gram (Parent & Tanis, 2011). 

More recently, the Post-School Achievement 
Through Higher Skills (PATHS) curriculum was 
developed to promote social cognitive career 
and self-determination outcomes and address 
the unique learning needs of young women with 
disabilities (Lindstrom et al., 2013). Lesson top- 
ics, derived from special education and career 
development literature, include self-determina- 
tion and strengths, gender identity, disability 
knowledge, and career exploration. Participants 
in PATHS were young women identified by 
school personnel as either eligible for special 
education services or at risk for dropping out of 
school (e.g. chronic absenteeism, credit defi- 
cient, unstable family situations). A pretest-post- 
test control group design with 110 young women 
found that that participation in this curriculum 
resulted in significant gains in autonomy and 
disability and gender-related knowledge (Doren 
et al., 2013). The Paths 2 the Future (P2F) cur- 
riculum used in the current study is an extension 
of the PATHS curriculum, with updated content 
and language to reflect current career- and col- 
lege-readiness standards. 

Overall, the reviewed literature highlights 
the potential benefits of gender-specific career 
and transition interventions; however, to date, 
there have been no randomized controlled 
studies targeting improved career outcomes 
for young women with disabilities receiving 
special education services. The current study 
builds on previous research to test the efficacy 
of a fully developed gender-specific career 
development intervention (i.e., P2F) using a 
randomized controlled trial. 


Social Cognitive Career 
Theory (SCCT) 


The P2F intervention utilizes SCCT as a frame- 
work for understanding the career development 
processes for young women with disabilities. 
SCCT is an integrated career development 
model that describes how people, their behavior, 
and environments influence each other to shape 
occupational and academic interests, choices, 
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and attainment of career goals (Lent, 2005). 
SCCT takes into account the influence of indi- 
vidual and contextual variables, such as race- 
ethnicity, gender, disability, and socioeconomic 
status, as well as self-efficacy and expectations 
for the future while at the same time accounting 
for the importance of significant learning experi- 
ences in shaping career behaviors and educa- 
tional outcomes. A key tenet of SCCT is that 
self-efficacy (beliefs about a person’s capabilities 
to achieve their goals) and outcome expecta- 
tions (imagined consequences of actions) have 
important effects on the formation of career 
interests, choices, and aspirations. Aspirations 
for the future, in turn, are a significant predictor 
of employment in adulthood (Ashby & Schoon, 
2010). Self-efficacy and outcome expectations 
predict career interests and work outcomes for 
students with academic learning challenges and 
intellectual disabilities (Nota et al., 2010; Ochs 
& Roessler, 2001), and previous interventions 
designed to enhance self-efficacy beliefs have 
been demonstrated to be successful for adoles- 
cents with disabilities (Lindstrom et al., 2019; 
Sheftel et al., 2014 ). Thus, we posit that self- 
efficacy and outcome expectations are crucial 
concepts for understanding and promoting 
career development outcomes among young 
women with disabilities. 


A key tenet of SCCT is that self- 
efficacy (beliefs about a persons 
capabilities to achieve their goals) 
and outcome expectations 
(imagined consequences of actions) 
have important effects on the 
formation of career interests, 
choices, and aspirations. 


Thus, we posit that self-efficacy and 
outcome expectations are crucial 
concepts for understanding and 
promoting career development 
outcomes among young women with 
disabilities. 


Study Purpose 


The purpose of this study was to determine 
whether the P2F curriculum produces a benefi- 


cial impact on career development outcomes for 
young women with disabilities relative to a 
counterfactual (e.g., typical high school career 
and transition services ). P2F is a fully devel- 
oped and pilot-tested curriculum designed to 
address the unique career development and tran- 
sition needs of young women with disabilities 
(Lindstrom et al., 2019). Developed and field- 
tested initially as PATHS (Lindstrom et al., 
2013), the P2F curriculum includes 75 lessons 
divided into four modules: (a) self-awareness, (b) 
disability knowledge, (c) gender identity, and (d) 
career and college readiness. The modules are 
designed to build individual skills and attributes 
and expose young women to a wide range of 
career options. Lessons are taught in sequence 
within a group or classroom setting in the course 
of the school day. (For more details on the P2F 
curriculum, see Lindstrom et al., 2019.) 

Guided by the SCCT framework and prior 
findings, we hypothesized that young women 
in the P2F intervention schools would exhibit 
greater rates of increased career development 
skills relative to those in the control schools. 
Hypotheses were tested using an intent-to- 
treat (ITT) approach. That is, effectiveness of 
a treatment relies on comparison of assigned 
groups regardless of dropout, changes in 
school protocol, variation in implementation 
or noncompliance of teachers, or any other 
unobserved factors that occur after random 
assignment (DeGarmo & Gewirtz, 2019). 
Thus, the randomization and ITT evaluation 
provide causal inference and unbiased esti- 
mates of intervention effects due to group 
assignment within real-world conditions. 


Method 
Study Procedures 


After receiving institutional review board 
approval, students were recruited from 26 pub- 
lic high schools in the northwestern region of 
the United States to participate in an efficacy 
trial of an intervention designed to improve 
career outcomes for young women with dis- 
abilities. Special education teachers and school 
counselors at each school were asked to iden- 
tify a sample of young women in their schools 
to participate in the study based on the follow- 
ing inclusion criteria: (a) identified as female; 
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(b) currently enrolled in Grades 9 through 12 in 
a participating high school; (c) determined eli- 
gible for special education services with a 
high-incidence disability, including learning 
disability, other health impairment (including 
attention deficit hyperactivity disorder), speech 
or language disability, and emotional disabil- 
ity; and (d) demonstrated fourth- to fifth-grade 
reading, writing, and language skills. Students 
with more significant cognitive or intellectual 
disabilities were not included. 

Schools were randomly assigned to inter- 
vention and control conditions using propen- 
sity score matching (PSM) procedures. PSM is 
recommended for group randomized trials of 
20 or fewer or correlations of .30 and above on 
matching factors (Diehr et al., 1995; Luellen 
et al., 2005; Murray, 1998). Over half of the 
matching factor correlations were above .30. 
Schools were matched on data extracted from 
the Department of Education on school size 
(number enrolled), percentage special educa- 
tion, percentage Black, percentage White, per- 
centage Hispanic, and percentage receiving 
free and reduced lunch. The propensity score 
(PS) was defined as the conditional probability 
of a subject being assigned to the treatment 
group given the observed covariates and where 
the log odds of assignment to treatment was 


inf =) =f, +B,School Size 


+B, %SPED+B,%EA 
+B, %AA +B, %HA +B,%FRL 


2 


where covariates were the size of the school and 
the percentage of students enrolled who were 
receiving special education, were white, were 
Black, were Hispanic, and were eligible for free 
or reduced lunch. The matching tolerance was 
set to .2, and paired scores were randomly allo- 
cated to either treatment or control. Participants 
in the control condition received “business-as- 
usual” career and transition services within their 
high school. A total of 156 participants from 13 
schools were randomized to the intervention 
condition, and 230 participants from 13 schools 
were randomized to the control condition. 
Parent consent and student assent were 
obtained prior to the study onset. Teachers 


were trained in the data collection instruments 
and procedures to be used during the study. 
Students in both control and intervention 
groups completed a web-based survey at four 
time points: (a) preintervention (T1), (b) mid- 
way through the academic year (for schools 
implementing a full-year schedule; T2), (c) 
postintervention (T3), and (d) 6-month follow- 
up (14). The P2F survey was administered 
individually or in groups. Each student was 
provided with a unique login code. Students 
took between 15 and 60 min to complete the 
survey. Research team members were avail- 
able for the entire duration to answer any ques- 
tions the students had about survey items or 
complex vocabulary. Teachers separately pro- 
vided information about each participating 
student’s diagnosis qualifying them to receive 
special education services as well as informa- 
tion on academic, family, health, work, and 
other barriers that often limit educational 
attainment and may influence postschool out- 
comes. 


P2F intervention. Prior to the start of the inter- 
vention and data collection activities in each 
school, all P2F instructors received one full 
day of professional development from the pro- 
gram developers and research team. This 
included an explanation of the P2F logic 
model, an overview of the content of the four 
program modules, and modeling of procedures 
for implementing the standard daily lesson 
plans. Once appropriate consent and assent 
processes were complete, intervention teach- 
ers began implementing the lessons and activi- 
ties as outlined in the P2F curriculum guide 
(Lindstrom & Post, 2015). Lessons included 
small- and large-group discussions, role-plays 
to practice key skills, guest speakers, and field 
trips to explore a wide range of college and 
career options. Young women attended the 
class daily (or on alternate days for those on 
block schedules) and received high school 
credit upon completing the 75 lessons. 

The research team used multiple measures to 
observe and document fidelity of implementa- 
tion for the P2F study. First, we collected infor- 
mation on student dosage through tracking 
student attendance and lesson completion at 
each school. Curriculum dosage was calculated 
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at a school level as the number of lessons deliv- 
ered divided by the total number of curricular 
lessons (75). Second, we conducted a minimum 
of eight classroom observations in each of the 
intervention sites and completed a five-item 
teacher adherence observation protocol and six- 
item quality-of-instruction observation proto- 
col. Both the adherence and quality observation 
rubrics were completed by P2F research per- 
sonnel during the same classroom observation. 
Across the 13 intervention schools, lesson dos- 
age ranged from 49% to 100% (mean dosage = 
87.8%), with average observer-rated quality at 
79.7% and overall adherence to lesson compo- 
nents measured at 81.7%. 


Sample. The current study sample comprised 
366 participants (mean baseline age = 16.54 + 
1.12 years) who had baseline data available 
on all outcome variables. The study sample 
(N = 366) included 61.5% non-Hispanic 
White and 19.4% Latina youth, enrolled in 
Grades 9 (14.5%), 10 (33.9%), 11 (30.1%), 
and 12 (21.6%). Most participants were receiv- 
ing special education services under the catego- 
ries specific learning disability (55%) or other 
health impairment (14.8%). (See additional dis- 
ability demographics in Table 1.) Teachers were 
asked to indicate whether each student was 
experiencing any additional risk factors or bar- 
riers in five areas, including academics, family 
or living, work, at-risk behaviors, and health 
challenges. Most common_ teacher-reported 
barriers included difficult family circumstances 
(43.7%), mental health issues (42.6%), chronic 
absences (27.6%), no prior work or volunteer 
experience (25.1%), or being behind in com- 
pleting credits toward graduation (16.4%). 
According to teacher reports, 83 participants 
(22.7%) did not experience any barriers, but the 
majority of participants (53.6%) experienced 
more than one barrier. 


Measures. We developed a web-based survey 
that included a compilation of validated mea- 
sures. The survey was originally developed and 
validated through a pilot study (Doren et al., 
2013) and was augmented for the current assess- 
ment. The survey included between 107 and 128 
questions, depending on responses to items 
with branching. In addition to demographic 


Table |. Paths 2 the Future Student 
Demographics. 


Demographic n % 
Grade 
9 53 14.5 
10 124 33.9 
11 110 30.1 
12 79 21.6 
Missing 20 5:2. 
Hispanic or Latina 
Yes 7\ 19.4 
Missing 20 5.2 
Race-ethnicity 
White 225 61.5 
Black or African 17 4.46 
American 
American Indian or 16 44 
Alaskan Native 
Asian or Pacific Islander 8 2.2 
Other 44 12.0 
More than one race 49 13.4 
Missing 7 1.9 
IDEA classification 
Specific learning disability 202 55:2: 
Other health impairment 54 14.8 
Emotional disturbance 22 6.0 
Intellectual disability 22 6.0 
No disability specified 21 5.7 
Speech or language 14 3.8 
impairment 
Other 25 6.8 
Missing 6 1.6 


Note. N = 366. IDEA = Individuals With Disabilities 
Education Act. 


questions (14 items), the survey contained 
questions on future goals (four items), career 
and technical education course enrollment (one 
to three items), work experience (one to 20 
items), and 87 scale items. Eight previously 
validated scale measures were used to compute 
the focal career development score. 

Two scales were from the Arc’s Self-Deter- 
mination instrument (Wehmeyer & Kelchner, 
1995). The Autonomy scale included 14 items 
rated from | (not even if I had a chance) to 4 
(every time I have the chance). Sample items 
were “I work on schoolwork that will improve 
career chances,” “I keep my appointments and 
meetings,” and “I volunteer for things I am 
interested in.” Cronbach’s alpha for internal 
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consistency was .74, .85, .82, and .83 for T1, 
T2, T3, and T4, respectively. The Self-Real- 
ization scale included 15 items rated from 1 
(never agree) to 4 (always agree). Sample 
items were “It is better to be yourself than 
popular,” “I feel free to be angry at people I 
care for,” and “I don’t accept my own limita- 
tions” (reversed) (a = .78, .78, .80, and .79 
for T1 through T4, respectively). 

The Self-Advocacy scale was from the 
College Students With Disabilities Campus 
Climate Survey (Lombardi et al., 2011) and 
included five items rated from 1 (strongly dis- 
agree) to 5 (strongly agree). Sample items 
were “I know about my rights and responsi- 
bilities as a student with a disability” and “TI 
feel comfortable advocating for myself and 
my needs at this school” (a = .79, .82, .84, 
and .80 for T1 through T4, respectively). 

Two scales employed in the pilot study 
were adapted from McWhirter et al. (2000), 
the Vocational Skills Self-Efficacy scale and 
the Vocational Outcome Expectancy scale. 
Vocational Skills Self-Efficacy included 29 
items measuring confidence in job prepara- 
tion skills, time management, and goal set- 
ting rated from 1 (no confidence at all) to 5 
(complete confidence). Sample items were 
“Know what to wear for a job interview,” 
“State my general career interests,” and “State 
my educational goals” (a = .97, .97, .97, and 
.98 for Tl through T4, respectively). The 
Vocational Outcome Expectancy _ scale 
included six items rated from 1 (strongly dis- 
agree) to 4 (strongly agree). Sample items 
were “My career planning will lead to a sat- 
isfying career for me” and “My talents and 
skills will be used in my career/occupation” 
(a = .88., .92, .91, and .89 for Tl through T4, 
respectively). 

Two scales were from the Student Engage- 
ment Instrument (Appleton et al., 2006). 
Future Aspirations and Goals included 5 
items rated from | (strongly disagree) to 4 
(strongly agree). Sample items were “Going 
to school after high school is important” and 
“T plan to continue my education following 
high school” (a = .88, .91, .91, and .91 for 
Tl through T4, respectively). Peer Support 
for Learning included six items. Sample 
items were “Other students here like me the 


way I am,” “Students at my school are there 
for me when I need them,” and “Students 
here respect what I have to say” (a = .90, .91, 
.93, and .92 for T1 through T4, respectively). 

Finally, Disability and Gender Awareness 
(Doren et al., 2013) was a six-item scale rated 
from | (no confidence at all) to 5 (complete 
confidence). Sample items were “Identify dif- 
ferent types of disabilities,’ “Describe how 
being female can impact career choices,” and 
“Identify key characteristics of women lead- 
ers” (a = .79, .87, .89, and .91 for T1 through 
T4, respectively). 

After factor analyses (see next section), an 
overall career development construct was 
computed as a mean total score of the eight 
subscales drawn from above. Each compo- 
nent score was rescaled | to 5 and reflected to 
assess career skills and then averaged to com- 
pute the overall career development score 
(Cronbach’s alpha was .82, .82, .86, and .87 
for T1 through T4, respectively). 


Analytic Strategy 


Analyses were conducted in two stages. The 
first stage focused on the measurement of the 
criterion outcome construct of career develop- 
ment. In this preliminary stage, we conducted 
factor analysis, more specifically, principal 
components analysis (PCA) to evaluate the 
configural structure of the career development 
construct over time as measured by the con- 
stituent scale scores. The second phase of the 
analyses focused on evaluation of the P2F 
efficacy hypothesis, which involved multi- 
level latent growth curve analyses to address 
the repeated-measures structure of the data 
and to address the nonindependence of data 
observations due to the clustering of the data 
within students and schools. 


Career development construct factor analy- 
sis. We used several related approaches to 
build the career development construct score. 
First, we conducted exploratory factor analy- 
sis (EFA) of the scale score indicators. Spe- 
cifically, we used PCA across T1 through T4 
using the following recommended criteria for 
determining the factor structure of the focal 
construct outcome score (Henson & Roberts, 
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2006; Thompson & Daniel, 1996): (a) PCA 
was conducted with oblimin rotation allowing 
exploratory factors to be potentially corre- 
lated, (b) eigenvalues of unrotated factors 
were greater than 1, and (c) each factor 
accounted for greater than 5% of the total 
variance. 

Second, to best determine the optimal 
number of factors we conducted both parallel 
analysis (PA) and Velicer’s minimum average 
partial correlation analysis (MAP) (O’Connor, 
2000). The PA directly compares eigenvalues 
for the ranked ordered PCA loadings with 
eigenvalues generated from randomly ordered 
data (O’Connor, 2000). Criteria for meaning- 
ful factors are those with higher eigenvalues 
from those randomly generated (i.e., eigne- 
values greater than the random parallel 95% 
quantile are retained), and thus meaningful 
variance is obtained above chance. The ratio- 
nale underlying MAP is that the average of the 
squared partial correlations should be mini- 
mum when an adequate number of factors are 
partialed or extracted. 

Finally, for tests of measurement invari- 
ance, we conducted a replicability test recom- 
mended by Osborne and Fitzpatrick (2012) in 
which the standardized structure item load- 
ings from wave to wave are compared using a 
longitudinal panel design. Squared differ- 
ences less than .04 represent replication; dif- 
ferences greater in magnitude do not replicate. 


Multilevel P2F efficacy analyses. This study 
used a clustered randomized control trial to 
test the efficacy of the intervention (Murray, 
1998; Raudenbush, 1977). The clustering, or 
nesting of data, involved longitudinal data 
nested within students and students nested 
within schools. The P2F efficacy hypotheses 
were tested with structural equation modeling 
(SEM) specified as latent growth models 
(LGM) using Mplus 8.2 (Muthén & Muthén, 
1998-2018). We employed multilevel SEM 
(MLSEM) to estimate the two-level growth 
model at the within and between school levels 
(Heck & Thomas, 2015). The MLSEM growth 
model statistically addresses the nonindepen- 
dence of the clustered data and partitions vari- 
ance in participants’ growth trajectories across 
student- and school-level effects. 


LGM provide advantages for testing 
developmental changes in student outcomes 
within randomized trials (Brown et al., 2008) 
because they more reliably estimate within- 
and between-student differences in change 
compared to repeated-measures ANOVA. In 
LGM, a two-level growth model is a special 
case of an SEM measurement model, with 
repeated-measures observations specified as 
the measurement indicators of a latent vari- 
able intercept and latent growth rate or slope 
factor. The first step of LGM is to estimate a 
series of sequential unconditional models to 
determine the optimal pattern of growth. 
Using fixed chronometric time weights, we 
compared three models: a random intercept- 
only model (i.e., no mean growth but indi- 
vidual differences or variance in intercepts), a 
linear growth model (i.e., mean increases or 
decreases and variance in individual trajecto- 
ries), and finally, accelerated or quadratic 
growth to capture curvature in the mean tra- 
jectories. The intercept-only model is speci- 
fied by fixing repeated waves at 1, 1, 1, and 1 
for Tl through T4. The linear growth slope 
model adds a latent factor with loadings fixed 
at 0, 1, 2, and 3, representing variation in lin- 
ear trajectories. Finally, a quadratic model 
adds a latent factor allowing for acceleration 
in growth in trajectories, or nonlinear curva- 
ture, by fixing loadings at 0, 1, 4, and 9. 

As in SEM, the MLSEM model fit was 
evaluated using recommended fit indices 
(Byme, 2011; McDonald & Ho, 2002) of a 
chi-square minimization p value greater than 
.05, a comparative fit index (CFI) greater than 
.95, a chi-square ratio (y7/df) less than 2.0, 
and a root mean square error of approximation 
(RMSEA) less than .08. Data were also mod- 
eled using full-information maximum likeli- 
hood (FIML), which uses all available 
information from the observed data in han- 
dling missing data. FIML estimates were 
computed by maximizing the likelihood of a 
missing value based on observed values in the 
data. Any individual who had baseline data 
only and no follow-up data contributed noth- 
ing to the likelihood of estimates and was 
effectively excluded from change analyses 
(Brown et al., 2008). Compared to mean- 
imputation, listwise, or pairwise models, 
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FIML provides more statistically reliable 
standard errors. 


Within-school two-level growth Level |. The 
within-school growth model, or student-level 
growth model, for career development was 
specified as the following: 
Levell1DC;, = No; +11; Time, 


“2 
+1), Time’, +€ 


tii * © rij? 


where the Level | career development out- 
come variable is a function of the repeated- 
measures random intercepts (1,) for student 
i in school j and is a function of time as linear 
growth (n,) and time squared as accelerated 
or quadratic growth (n,) for student 7 in school 
j, plus the time-specific errors of measure- 
ment nested within individuals and schools 
(€,,). Thus Level 1 includes the time-varying 
repeated measures of the outcome score with 
linear time weighted 0, 1, 2, and 3 for the 
respective outcome scores and time squared 
weighted 0, 1, 4, 9. 


Within-school growth Level 2. The estimated 
time-varying random intercepts and growth 
trajectories for each student are then summa- 
rized for students representing the within- 
school but “between-student” effects specified 
as the following: 


Level2 CD Random 
Intercept), =O); + Soy 
Level2 CD Linear 
Growth n,,, = 0, + ¥,;Grade, 
+¥,, Ethnicity, + y,,SES, +4), 

Level2 CD Accelerated Growthn,,, = 

a,, +7, Grade, + y,,, Ethnicity, 

+75, SES, +0 5p 


where 1),, is the summary of individual differ- 
ences in student intercepts at Wave | for the 
full sample, Ni is the predicted growth rate 
for the sample, o.,, is the mean adjusted growth 
rate, and y,,, through y,,, are the student-level 
covariate effects of girls’ grade, ethnicity, and 


family socioeconomic status, respectively, as 
predictors of growth in career development. 
The residual errors for initial status, linear and 
quadratic growth (Coy Cup and € >) are each 
assumed to be normally distributed. 


Between-school model Level 3. The between- 
school growth model estimates the proportion 
of variation in latent intercepts and latent 
growth factors that is attributable to school- 
level variables. The latent growth model is 
identical to the specification of time weights 
at the within-school student level, and the 
growth indicators are latent estimates based 
on the summary of the estimated observed 
data for students in the within model. All 
school-level covariates predicting growth are 
entered in the between-school model. More 
importantly, because randomization occurred 
at the school level, the ITT P2F effect is esti- 
mated in the between model as 


Level3CD Interceptn, ; =a), +o; 
Level3CD Linear Growth n, , = a, ; +B, , P2F; 
+B, SchoolSize, +B,;,;%FRL, +6), 
Level3CD Accelerated Growth n, , = a,, 
+B,,P2F, +B, ,SchoolSize, 
+B,, %FRL, +¢,, 


where By is the ITT effect of the P2F inter- 
vention on growth in girls’ career develop- 
ment scores, controlling for the effects of 
school-level variables, such as school size 
(B,,) and percentage of free and reduced 
lunches served (B;)). The P2F efficacy 
hypothesis in the form of an MLSEM growth 
model is summarized graphically as a two- 
level within-school portion and a two-level 
between-school portion of variance, shown 
in Figure |. For clarity’s sake, model specifi- 
cation is illustrated with linear growth 
parameters only. The observed data or mani- 
fest variables are represented as squares. All 
unobserved latent variables are represented 
by circles, such as Level 2 and Level 3 
growth factors and the between-school latent 
estimates of career development repeated- 
measurement indicators. 
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Initial Status 
Intercept 


Between 
Initial Status 
Intercept 


Within School Model 


Figure |. Multilevel structural equation modeling analytic specification for test of the Paths 2 the 
Future efficacy hypothesis predicting growth in young women’s career development outcomes. 


Results 


Career Development Construct 


The means and standard deviations for the 
career development indicator scale scores are 
presented in Table 2 by group condition 
across time. In the first step of the analyses, 
we subjected each of the candidate scale 
score outcomes to a PCA, followed by the PA 
and the MAP focusing on the optimal number 


of factor dimensions. Results of the PCA, PA, 
and MAP each indicated that a single-factor 
solution was the optimal configural structure 
for the career development construct. Load- 
ings for the obliquely rotated PCA solution 
are shown in Table 3 for T1 through T4. 

For Tl, T3, and T4, the PCA obtained a 
single-factor solution and no second factor 
was extracted meeting criteria. At T2, the 
assessment wave with partial data, the PCA 
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Table 4. Means and Standard Deviations for Career Development Indicators by Group Condition. 
Control condition P2F intervention 

Survey administration n M SD n M SD 
Time | 213 3.43 0.45 153 3.38 0.57 
Time 2 4| 3.44 0.57 46 3.53 0.48 
Time 3 192 3.44 0.53 136 3.62 0.57 
Time 4 165 3.43 0.56 102 3.59 0.56 


Note. Scale score indicators were rescaled | to 5 and averaged to compute career development construct score. 


obtained two factors; however, only the first 
factor represented meaningful variance 
according to the PA test, in which the first 
factor eigenvalue was the only factor higher 
than randomly ordered variables. That is, the 
second factor was not meaningful variance 
above chance. Similarly, Velicer’s MAP test 
indicated that only the first factor had ade- 
quately minimized squared partial correla- 
tions. That is, when comparing component 
variance relative to error variance, the partial 
correlations will be small; as more compo- 
nent variability is removed or partialed out, 
the error variance will be larger and the par- 
tial correlations rapidly increase. 

In the last step of the factor analyses, we 
examined the replicability of the factor struc- 
ture over time (i.e., configural invariance). 
Results of the Osborne and Fitzpatrick (2012) 
replicability indicated that none of the squared 
PCA loadings differed among study waves 
across time. That is, when using the loadings 
shown in Table 3, comparing the squared abso- 
lute difference for all possible paired compari- 
sons revealed that none of the differences were 
greater than .04. 


P2F Efficacy Hypothesis 


The first step of the efficacy hypothesis was to 
determine the optimal pattern of change in the 
focal career development construct for the 
sample. We first inspected the means and 
standard deviations by group condition pre- 
sented in Table 4. The data in Table 4 indi- 
cated that young women in the control schools 
exhibited a rather flat trajectory over time, 
whereas those in the P2F intervention schools 
exhibited increases in means across time. 


We next estimated an unconditional multi- 
level growth curve analysis of the repeated- 
measure career development scores (see 
Figure 1). We examined a random intercept 
model only, a model adding in a linear slope 
factor, and finally, a model adding a nonlinear 
quadratic factor. The best-fitting model was 
the latter model including both a linear slope 
and a quadratic slope. Mean factor scores are 
estimated at the between-school level, and 
individual student variation in slope trajecto- 
ries are estimated at the within-school level. 
The means and variances for the uncondi- 
tional model are presented in Table 5. For the 
between-school means, the initial status inter- 
cept was significantly different from zero 
(M = 3.408, p < .001) as were the mean lin- 
ear growth rate (M = .107, p < .01) and the 
mean quadratic growth (M = —.027, p < .05). 
The negative mean for the accelerated growth 
factor indicated the shape of the curvature was 
a downward-cupped shape meaning that on 
average, the sample increased in the first year 
and then leveled off or decelerated by the sec- 
ond school year. The unconditional quadratic 
growth model obtained excellent fit to the 
data, y°(1) = .274, p = .601, CFI = 1.00, 
RMSEA = .00; ¥7/df = .45. Inspection of the 
within-school variances showed that there 
were individual differences in initial-status 
career development (oc? = 2.15, p < .001); 
however, the variance for the linear and qua- 
dratic factors were not different from zero, 
meaning that students grew in similar linear 
trajectories and similar patterns of nonlinear 
trajectories. To illustrate the dependent crite- 
rion factors for the efficacy hypothesis, the 
observed and estimated means from the multi- 
level SEM growth curve model are plotted in 
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Table 5. Unconditional Latent Factor Means and Variances for Multilevel Growth in Career 


Development. 


Career development 


Between school Within school 
Model M SE Variance SE 
Random intercept 3.408*** .026 215Re .050 
Linear growth slope 0.107** .040 059 .039 
Quadratic growth slope -0.027* 012 .004 .006 
Note. Model fit, x2(1) = .274, p = .601, comparative fit index = 1.00, root mean square error of approximation = .00; 


yidf = 5. 
*p < 05. **p <.01. **p < 001. 


———o——— Observed Means Control Schools 
---4--- Estimated Means Control Schools 
———o—— Observed Means P2F ITT Schools 
=== 4- == Estimated Means P2F ITT Schools 
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Figure 2. Observed mean trajectories and estimated mean trajectories by group condition for best- 
fitting multilevel structural equation modeling latent growth model including a random intercept, linear 
growth, and quadratic or accelerated growth factors, y7(1) = .274, p = .601, comparative fit index = 1.00, 
root mean square error of approximation = .00; y2/df = .45. 


Figure 2 for the control and P2F intervention 
schools. The plot indicates that a linear growth 
trend adequately fits the control schools, 
whereas a linear increase with a negative qua- 
dratic, or decelerated leveling off, best describes 
the observed and fitted data for the intervention 
schools. Because there was no significant varia- 
tion in linear or quadratic growth, the variance 
for the between-school growth factors are fixed 


at zero (Heck & Thomas, 2015), whereas the 
intercept factor variance is freely estimated. 
In the final step of the analyses, we tested the 
ITT efficacy hypothesis specifying the effect of 
the randomized group contrast at the between- 
school level controlling for school characteris- 
tics of size, rurality, percentage free and reduced 
lunch, and type of academic term implementa- 
tion of the curriculum (e.g., trimester, semester, 
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Within School Model 
Growth Career Development 
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Initial Status 
Intercept 


Accelerated 
Growth 
Slope 


Between 
Initial Status 
Intercept 


Between 
Linear 
Growth 

Slope 


Between 
Accelerated 
Growth 
Slope 


Between School Model 
Growth Career Development 


Figure 3. Career development multilevel latent growth model test of Paths 2 the Future intent-to- 
treat efficacy hypothesis. Estimates are unstandardized coefficients; y7(21) = .004, p = 1.00, comparative 
fit index = 1.00, root mean square error of approximation = .00; 7/df = .00. 

tp <.10. *p <.05. **p < 01. ***p < 001. 


full-year block schedule). Student-level covari- _ unstandardized regression coefficients. As sug- 
ates were entered at the within-school level. gested in the mean plots, the efficacy hypothesis 
Results of the conditional MLSEM growth — was supported. Controlling for student and 
model are presented in Figure 3 inthe form of —— school characteristics, we found that young 
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women in the P2F intervention schools grew in 
career development skills at a greater rate rela- 
tive to those in the control schools (8 = .18, p 
< .01). The P2F ITT contrast was also associ- 
ated with greater deceleration relative to the 
controls (8 = —.05, p < .05). Among the 
school-level covariates in the between-school 
model, free and reduced lunches were associ- 
ated with higher levels of career development 
at preintervention baseline Tl (B = .37, p < 
.01). Because the growth factor variances are 
fixed at zero, covariances (correlations) of the 
latent factors are not estimated at the between- 
school level. 

Among the  student-level predictors of 
growth, the teacher-rated student risk index was 
associated with lower levels of career develop- 
ment skill at baseline (y = —.05, p < .001), 
and socioeconomic status was marginally asso- 
ciated with higher levels of career development 
skill (y = .06, p < .07). Among the growth fac- 
tors, initial status was negatively associated with 
linear growth and positively associated with 
accelerated growth, meaning those students 
lower in career development at baseline were 
more likely to increase, and those initially higher 
grew at a slower rate. In other words, findings 
suggest that participants with higher levels of 
teacher-reported barriers had lower career devel- 
opment scores at baseline, and those with lower 
baseline scores showed more improvement over 
time. Overall, the MLSEM model obtained 
excellent fit to the data, y7(21) = .004, p = 1.00, 
CFI = 1.00, RMSEA = .00; 7/df = .00. 

In sum, data reduction of the primary 
career development scale scores indicated that 
scale scores adequately defined a single-fac- 
tor solution for the career development con- 
struct. Furthermore, the factor structure was 
consistent over time. The ITT analyses sup- 
ported effectiveness of the P2F intervention 
for young women with identified disabilities. 
Effect size of the P2F intervention was calcu- 
lated using Cohen’s d and from the Y standard- 
ized effects in a single-level growth model 
to obtain student-level effect size. Using the 
means in Table 5, Cohen’s d was .33 at T3 and 
.29 at T4, indicating a moderate effect. Cohen 
(1988) characterizes a d of .2 as small, .5 as 
medium, and .8 and above as large. The stu- 
dent-level effect size from a Y standardized 


solution was a medium effect of the interven- 
tion predicting growth (effect size [ES] = 
.56). The effect is larger because of the 
increased reliability in estimating change over 
time as compared with cross-sectional mean 
comparisons. Finally, the between-school ES 
was large (ES = 1.7) according to the Y stan- 
dardized solution from Figure 3. The between- 
school effect is very large due to the large 
number of schools (26) relative to the 366 
young women at the student level. 


Discussion 


This study tested the efficacy of the P2F 
career development curriculum, designed to 
address the unique barriers faced by young 
women with disabilities in transition from 
high school to college and careers. Data 
reduction of the primary scale scores ade- 
quately defined a single-factor solution for the 
multifaceted career development construct 
that included self-efficacy, outcome expecta- 
tions, and disability and gender awareness 
scales, among others. The factor structure was 
also consistent over time, suggesting that the 
construct mean change between groups repre- 
sented reliable and valid differences attributed 
to the P2F intervention. In a unique contribu- 
tion to the literature, we tested the efficacy 
hypothesis through an ITT analysis of this 
gender-specific intervention. Controlling for 
student and school characteristics, we found 
that young women in the P2F intervention 
schools grew in career development skills at a 
greater rate relative to those in the control 
schools. The student-level effect size was 
medium, and the between-school effect size 
was large. Together, the data suggest that P2F 
was effective in promoting individual student 
differences in career development outcomes 
for young women with disabilities who par- 
ticipated in this study. 


Together, the data suggest that P2F 
was effective in promoting 
individual student differences in 
career development outcomes for 
young women with disabilities who 
participated in this study. 
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Young women with disabilities face mul- 
tiple complex barriers during the transition 
from high school to adulthood, and their 
career options and opportunities are often 
influenced by the confluence of gender ste- 
reotypes and disability barriers (Lindstrom 
et al., 2012; Hogansen et al., 2008). Despite 
the well-documented gaps in career expecta- 
tions and outcomes for this population, P2F 
is one of very few comprehensive curricula 
designed and tested to provide specific 
instruction in career development concepts 
for young women with disabilities, including 
lessons on topics of self-awareness, disabil- 
ity knowledge, gender identity, and career 
and college readiness. 

Although other career and transition pro- 
grams include instruction on self-determina- 
tion or disability awareness (e.g., Sheftel et al., 
2014; Wehmeyer et al., 2009), P2F is unique in 
that participants are introduced to both disabil- 
ity and gender content to inform career aware- 
ness and expand career exploration. After 
exploring their personal strengths, young 
women in P2F participate in a set of lessons 
and activities focused on disability knowledge, 
educational rights and responsibilities, self- 
disclosure, and accommodations. With a 
clearer understanding of disability identity, 
participants are then exposed to lessons explor- 
ing gender roles, expectations, and stereotyp- 
ing and learn about barriers and opportunities 
for women in the workplace. We believe that 
the unique combination of disability and gen- 
der lessons in the P2F curriculum influenced 
increases in disability and gender awareness 
for intervention participants and also impacted 
the overall career development outcomes in 
this study. Our findings regarding improve- 
ments in gender and disability awareness also 
confirm and extend an earlier quasiexperimen- 
tal study of the PATHS curriculum, a precursor 
to P2F (Doren et al., 2013). 

Utilizing SCCT as a guiding theoretical 
framework (Lent, 2005), we were also able to 
examine important individual variables, such 
as self-efficacy and outcome expectations, 
across our sample. Previous studies have found 
that young women with disabilities in high 
school often have low self-esteem and feel less 
confident than their male peers about entering 


high-wage or nontraditional occupations (Ferri 
& Connor; 2010; Lindstrom et al., 2012). Ina 
survey of adolescents with disabilities and their 
families, Powers and colleagues (2008) found 
that females were more likely than males to 
indicate that people expect less of them based 
on their disability, and both females and par- 
ents of daughters reported constrained expecta- 
tions due to gender. Despite the persistent and 
pervasive low expectations for this population, 
young women participating in P2F were 
exposed to new content that led to important 
gains in outcome expectations and vocational 
skills self-efficacy, essential constructs that 
have been shown to influence the formation of 
career aspirations and choices (Lent, 2005). 

Finally, we were also interested in under- 
standing the influence of a gender-specific cur- 
riculum (P2F) on individual attributes and 
skills, such as autonomy, self-realization, and 
self-advocacy. Previous studies have docu- 
mented the importance of these core compo- 
nents of self-determination in promoting the 
school- and career-related adjustment of ado- 
lescent girls with disabilities or other risk fac- 
tors (Doren et al., 2013). Researchers who 
have examined the impact of mixed-gender 
career development programs have found that 
targeted career development instruction can 
lead to significant increases in self-realization 
and autonomy for both high school students and 
adults with disabilities (e.g., Sheftel et al., 2014; 
Shogren et al., 2016); however, these programs 
did not also attend to gender or disability. Our 
findings documenting growth in career devel- 
opment skills for young women participating in 
the gender-specific P2F curriculum extend the 
literature and add to the limited body of research 
that considers both gender and disability in the 
context of career development. 


Limitations 


There following limitations should be taken into 
consideration when interpreting the results. First, 
teachers and school counselors in the interven- 
tion schools referred young women to participate 
in the study, which may have introduced some 
selection bias. Second, teachers and data collec- 
tors were not blind to their assignment to inter- 
vention versus control condition, which could 
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have resulted in potential performance bias. 
However, research team members were present 
at the time of all survey data collection and con- 
ducted fidelity observations on random occa- 
sions while the intervention was being delivered. 
Participants in the intervention condition were 
also aware of their assignment, and given our 
study design, we were unable to control for such 
placebo effects. Finally, although our outcome 
variables were assessed using self-report mea- 
sures, our use of individual web-based surveys as 
opposed to face-to-face interviews should have 
minimized self-report bias. 

Notwithstanding these limitations, this 
study adds to the sparse body of literature 
focused on career development for young 
women with disabilities enrolled in high 
school. Although previous studies have used 
quasiexperimental designs or mixed-gender 
career interventions, this was the first random- 
ized controlled trial testing the efficacy of a 
gender-specific career development curricu- 
lum for young women with disabilities. More- 
over, the school randomized clustered trial 
involved PSM of school characteristics prior 
to randomization. Therefore, in addition to 
randomizing schools, matching procedures 
further reduced threats to group comparisons 
and internal validity of the study design. In 
addition, we used latent growth models as 
described earlier. Growth models have become 
more standard in assessing multiple-wave data 
because they more reliably assess change than 
do pretest-posttest change models only (Dun- 
can et al., 2006; Singer & Willett, 2003). An 
advantage of the multilevel structure is that 
individual student differences in trajectories is 
modeled as variance predicted, not just group 
means of randomly assigned conditions. Fur- 
thermore, the trajectories account for nonlinear- 
ity in the data, providing greater reliability than 
linear assumptions alone. Brown et al. (2008), 
among others, have shown the advantages of 
this approach over repeated-measures ANOVA 
models in testing the efficacy of ITT designs. 


Implications for Research and 
Practice 


The findings from this efficacy trial have 
implications for both research and practice. 


Based on the results of this study, we found 
that participation in a gender-specific career 
development curriculum can have important 
and unique benefits for young women with 
disabilities. Given that young women with 
disabilities often have limited career aspira- 
tions and restricted career opportunities, we 
believe high schools special education pro- 
grams should focus on introducing content 
related to both disability and gender to 
encourage personal development and build 
awareness of a broad range of career oppor- 
tunities. We also recommend that teachers 
develop opportunities to reflect the unique 
needs of young women with disabilities in 
coeducational career and transition services, 
focusing on developing self-advocacy and 
self-awareness and increasing vocational 
skill self-efficacy. 

In terms of future research, it is crucial to 
continue to identify and explore gender dif- 
ferences in school experiences, transition 
planning, and postschool career outcomes for 
young women with disabilities receiving spe- 
cial education services. More specifically, 
single-sex career development interventions 
should be examined with more diverse sam- 
ples from various regions of the country, tak- 
ing into account race-ethnicity, language, 
and disability status, to understand potential 
differential impacts with a variety of par- 
ticipants and settings. Future studies could 
also employ a comparative effectiveness 
design to determine whether a gender-spe- 
cific career development curriculum is more 
effective than nonspecified general career 
development. In addition, we recommend 
that future studies in this area collect follow- 
up data either 1 or 2 years postintervention in 
order to document longer-term changes in 
career development constructs as well as 
enrollment in postsecondary education and 
engagement in employment. Finally, devel- 
opment and examination of school-based 
career development interventions for young 
men with disabilities are needed to help ele- 
vate the postsecondary career outcomes for 
all youth with disabilities. 

The P2F study extends previous career 
and transition research by testing the effi- 
cacy of a gender-specific career curriculum 
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for young women with high-incidence dis- 
abilities, such as learning disabilities and 
other health impairments, using a randomized 
trial. Unlike previous studies, P2F used a 
multidimensional career development con- 
struct that combined key indicators, includ- 
ing self-determination, disability and gender 
awareness, vocational skills self-efficacy, 
vocational outcome expectations, and student 
engagement. We found significant overall 
gains in career development over time, con- 
trolling for both student and school character- 
istics. Controlling for individual barriers, 
girls in P2F demonstrated greater gains rela- 
tive to the control group counterparts. Thus, 
we found that young women in this study 
with barriers that often impede high school 
completion or postschool success benefited 
significantly from participation in P2F. 
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